Theory of Fano Resonances in Graphene: The Kondo effect probed by STM 
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We consider the theory of Kondo effect and Fano factor energy dependence for magnetic impurity 
(Co) on graphene. We have performed a first principles calculation and find that the two dimensional 
E\ representation made of dxz , dyz orbitals is likely to be responsible for the hybridization and 
ultimately Kondo screening for cobalt on graphene. There are few high symmetry sites where 
magnetic impurity atom can be adsorbed. For the case of Co atom in the middle of hexagon of 
carbon lattice we find anomalously large Fano g-factor, 5 ~ 80 and strongly suppressed coupling to 
conduction band. This anomaly is a striking example of quantum mechanical interference related 
to the Berry phase inherent to graphene band structure. 

PACS numbers: 



Scanning tunneling microscopy (STM) allows us to 
probe the electronic properties of conducting materials 
with atomic scale spacial resolution. This experimental 
technique is particularly well suited to study electron- 
correlation phenomena. One of the most famous cor- 
relation phenomena is the Kondo effect arising from a 
localized magnetic moment being screened by the con- 
duction electrons^. It results in a sharp Abrikosov-Suhl 
resonance in local density of states (LDOS) of the impu- 
rity at the Fermi level and below a characteristic Kondo 
temperature, Tk- While the Kondo effect is well under- 
stood for impurities in bulk materials and simple model 
systems, STM has been substantial for revealing the in- 
tricacies of Kondo effect at conventional metal surfaces. 
One well known example is the Kondo effect caused by 
Co ad-atoms on a Cu (111) surface^: It shows that the 
rich electronic structure of three-dimensional metals like 
Cu, in general, makes impurity effects at their surfaces'^^^ 
depending strongly on atomistic details and requires to 
understand interaction mechanisms in detail. 

Graphene - a monolayer of carbon atoms arranged 
in a honeycomb lattice - is the first truly two dimen- 
sional material^ and provides a two-dimensional electron 
gas with distinct and highly symmetric low energy elec- 
tronic structure: At two non-equivalent corners of the 
Brillouin zone, K and K', the linearly dispersing valence 
and conduction band touch forming a conical point and 
leading to the Berry phase tt^'^. Thus, electronic ex- 
citations in graphene resemble massless Dirac fermions 
with the speed of light being replaced by the Fermi 
velocity Vi k, c/300. Therefore, graphene provides an 
important model system for understanding quantum ef- 
fects in reduced dimensions and in presence of an "ultra- 
relativistic" conduction electron bath. 

A theoretical study showed that even in undoped 
graphene the Kondo effect can exist above a certain crit- 
ical coupling despite the linearly vanishing density of 
states^^^"'^, a situation very similar to magnetic impurities 



in the pseudogap phase of high Tc superconductors^' 
Moreover, back-gating^ as well as chemical doping-'^'^'^^ 
allows one to control the chemical potential in graphene 
and to tune Kondo physics and electron tunneling in this 
way. 

In this paper, we address how the Kondo effect mani- 
fests in STM experiments on graphene and why Fano res- 
onances in the STM spectra can depend unusually strong 
on the chemical potential as well as the real space posi- 
tion of the impurity. To this end, we firstly consider the 
single impurity Anderson model with graphene provid- 
ing the host electronic structure and compare to simple 
model of usual metal surface. With this background, 
we turn to a more realistic ab-initio based description of 
magnetic impurities on graphene and discuss the case of 
Co ad-atoms as in the recent experiment by Manoharan 
et al.^'^. By comparison to Co on Cu (111), an extensively 
studied system"^"^ posessing also hexagonal symmetry of 
the surface, demonstrate the particular importance of im- 
purity induced resonances in graphene. Furthermore, we 
analyze impurities being bound to different sites of the 
graphene lattice and show that there one can expected 
a strong adsorption site dependence of Fano factors in 
STM experiments. 

Model for electron tunneling close to impurities. The 
TT-band the tight-binding Hamiltonian of graphene reads 



as 
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(1) 



where and hi are the Fermi operators of electrons in the 
carbon pz orbital of sublattice atoms A and B in the cell 
at Ri , respectively. The sum includes all pairs of nearest- 
neighbor carbon atoms and t k, 2.7 eV is the hopping 
parameter. With the Fourier transformed operators Ofc 
(6fc), defined by a; = /^^ j^^e^^^'^au and bi analogously. 
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the Hamiltonian can be rewritten as 







Of 
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— -^^iRk^k with *(fc) 
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Hk is the /c-dependent 2x2 matrix 





r(fc) 



(2) 
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with C(fc) = -tE^=ie'''^^^~^^\ bj {j = 1,2,3) are the 
vectors connecting neighboring atoms"'^^, and f^s is the 
area of the Brillouin zone. An impurity contributing a 
locahzed orbital, J?i,„p = e^^^p^^ d^^d^ + Un^ni, with 
Fermi operator d, energy esmp and on-site Coulomb repul- 
sion U is considered. Its hybridizations with the graphene 
bands is described by, V = J2k a "^1 cr^^da -I- h.c. This 
problem has been extensively discussed for normal metals 
and is usually called "Anderson impurity model" ^ . 

In this framework, the connection between a tip and a 
sample in the STM experiment can be expressed by the 
transfer Hamiltonian 



M = Y,{MdtdU 



-KH.c.)+V(*L^4ti.+H.c.) (4) 



describing tunneling of electrons from and to the STM 
tip with the tunneling matrix elements Mdt and M^t and 
the Fermi operators (ij.) for electrons in the STM tip. 

The Fano g- factor in the STM AI /dV spectra can be 
understood in terms of this model by using the equation 
of motion approach from Ref. 18. One finds 



A 
B 



with 



and 



A = Mdt + Y.MktVkP 



E-p — ty 



(5) 



(6) 



(7) 



where P is the principle value symbol. 

To obtain qualitative insights we proceed by simplify- 
ing these expressions: For a conduction electron state fc, 
denote its probability density integrated about an atomic 



with the 



local conduction electron Green function 
— As argued in Ref 5, this sim- 



-efc-iO+ ■ 

pie model has proved successful to describe Fano factors 
for CoCu„ clusters on Cu(lll) and will be used here to 
understand Fano resonances in graphene and why they 
are different to metals like Cu. 

Density of states effects on the Fano factor. Usually, 
the magnetic orbitals of the ad-atom are strongly local- 
ized resulting in \Mdt\ <C |Af | and consequently 



g « Re G(£;F)/Im G{Ey) 



(11) 



In a metal with bandwidth D and constant density of 
states (DOS) in the vicinity of the impurity, Im G(e) ~ 

-1x120 if -D < e < D, wc obtain q In 



D+Ef 
D-Ef 



2 Bp 



. Hence, \q\ < 1 and the Kondo effect manifests in 
STM on normal metals as anti-resonance close to Ep as 
long as \Mcit\ <C |Af |. This is very different for graphene: 
The graphene DOS is Nl{E) = B • Q{D - \E\) result- 



ing GiE) = ^ In 



E^ 



iTTN^iE) and 



2 sign (Ef) 



In 



Ef 
D 



(12) 



This results follows directly from linearity of N{E) and 
the Kramers-Kronig relations. As D « 6eV and usually 
Ep < 0.5 eV, the q-factor can be g ^ 1 and the Kondo 
effect may manifest in STM as resonance instead of an 
anti- resonance even for \Md,t\ This is in con- 

trast to a normal metal, where predominant tunneling 
into the conduction electron states results in a Kondo- 
antiresonance in STM. Moreover, Eq. (12) demonstrates 
that the q-factor in graphene can be expected to depend 
strongly on the chemical potential. 

Energy dependence of the asymmetry factor. Any im- 
purity being coupled to graphene leads to characteristic 
resonances in the local density of states in the vicinity 
of the impurity^^ and may consequently alter the Green 
functions to be inserted into Eq. (11). To understand the 
role of resonances in the local electronic structure for the 
q-factor, we illustrate the situation of a realistic impurity 
by comparing the experimentally important cases of Co 
on graphene and Co on Cu (111). 

For a realistic description of these systems wc per- 
formed density functional calculations within the gen- 
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spherecenteredattheCo-atomby l-^fcp. Then, ifthetip cralized gradient approximation (GGA) on 6 x 6 



is directly above the Co impurity, one can approximate 
MktVdk = |*fcpMX^ (8) 

with MV independent of k. Thus, Eq. (6) and Eq. (7) 
yield 



and 



A = Mdt + MVRe G(Ef) 



B = MVlm G{Ef) 



(9) 



(10) 



graphene supercells containing one Co ad-atom as 
well as on Cu(lll) slabs containing 5 Cu layers and 
one Co ad-atom. The Vienna Ab Initio Simulation 
Package (VASP)^^ with the projector augmented wave 
(PAW)^^'^"^ basis sets has been used for solving the re- 
sulting Kohn-Sham equations. In this way we obtained 
relaxed structures for both systems. In particular, we 
found that GGA predicts Co to sit above the middle of 
a hexagon on graphene. 

To estimate the Fano g-factors we extracted the or- 
bitally resolved Green functions at the impurity site using 
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FIG. 1; Hybridization functions A (a,c) and calculated asym- 
metry factors |q| for Co on grapliene (upper panel) and Co 
on Cu (lower panel). For the hybridization functions Re A is 
plotted as solid line, ImA dashed. 



atomic orbitals naturally included in the PAW basis sets: 
The projectors {di\ipn'k) of orbitals \di) localized at the 
impurity atoms onto the Bloch eigenstates of the Kohn- 
Sham problem j-fAnk) are available when using PAW as 
implemented in the VASP and these give the local Green 
functions according to 
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e + iS - e„k 



(13) 



So, we employ here the same representation of localized 
orbitals as used within the LDA+U-schcme implemented 
in the VASP-code itself or as discussed in the context of 
LDA-hDMFT in Ref. 24. 

The local Green functions at the impurity sites as de- 
fined in Eq. (13) are 5x5 matrices which can used to 
obtain the hybridization function A(e) of the impurity: 



G-\e) = e + tS-A{e). 



(14) 



Hence, A(e) are also 5x5 matrices describing hybridiza- 
tion of 5 d-electrons of Co. In the particular case of 
Co on Cu(lll) and graphenc. which arc both hexago- 
nal surfaces, these matrices are diagonal and decompose 
into degenerate blocks of two 2-dimensional and one 1- 
dimensional representations, transforming under the ro- 
tation group Cgi, as Ei, E2, and Ai. These components 
of the hybridization function are depicted in Fig. 1 a 
and c. At energies close to the Fermi level of graphene 
all graphene states are in the vicinity of the two Dirac 
points. These states transform under Cgt, according to 
El and E2. Hence, the hybridization of the Ai impu- 
rity orbital to the graphene bands is strongly suppressed. 
Moreover, the crystal field splitting appears to be such 
that the Ei orbitals {dxz and dyz) are highest in energy 
by approx 0.5-leV as compared to the other d-orbitals. 
So, the El orbitals are expected to determine the q-factor 
in STM experiments probing the Kondo effect of Co on 
graphene. 



For Co on graphene Fig. 1 shows that |ImA(e)| ^ 
Re A(e)| in the vicinity of e = 0, which is the Fermi level 
for undoped graphene. This is very different from the 
case of Cu, where |Im A(e)| and |Re A(e)| are mainly on 
the same order. Using 



e + i5 - Ek 



in combination with Eq. 
at 



(15) 

(8) and \Mdt\ < \M\ we arrive 



g « Re A(£;F)/Im A(£:f). 



(16) 



Within this approximation the projectors and eigencner- 
gies obtained from DFT allow for an ab-initio prediction 
of g-factors. The Fig. 1 b and d show the g-factors pre- 
dicted for channels of different Ce rotational symmetry 
as calculated for Co on graphene and Cu, respectively, 
as function of the resonant energy E. E = corresponds 
here to the Fermi level of the undoped system. The cal- 
culated g-factors for Co on Cu(lll) are typically on the 
order of g ^ 1 without pronounced energy dependence. 
This is in contrast to graphene, where g > 1 in a wide 
energy rage and q is strongly energy dependent. So, g is 
expected to be strongly sensitive to local changes in the 
chemical potential of graphene, which can be caused by 
gate voltages, chemical doping or substrate effects. We 
also point out that this discussion goes well beyond a 
simple tight binding and linearized dispersion analysis. 
Presence of defect can substantially change local bands 
and DOS. 

Site dependence of hybridization matrix elements. For 
Co on graphene, we saw that graphcne's Fermi surface 
being made up by states transforming as Ei and E2 under 
Cg lead to particular Co orbitals being decoupled from 
the graphene bands. This special symmetry of graphcne's 
Fermi surface makes the g-factors seen in Kondo res- 
onances in STS particularly dependent of the precise 
atomic arrangement of the magnetic impurity. This can 
be illustrated by the site dependence of g-factors for an 
Anderson impurity sitting , on top of C and on a bridge 
site in the middle of a hexagon, respectively. 

At each adsorption site, a spherically symmetric s-wave 
impurity can be modelled by equal hopping matrix ele- 
ment Vi to all adjacent sites. For such an Anderson im- 
purity on top of a carbon atom or at a bridge site we 
obtain 



and Vh = 



(17) 



respectively, by translating that nearest neighbor hop- 
ping into the matrix formalism of Eq. (3) and perform- 
ing the Fourier transformation. Combining this with Eq. 
(7) results in B - 0{Ef) for Ep 0. This situation 
corresponds to Eq. 12 with g being enhanced as E-p — > 0. 

For the impurity in middle of the hexagon the Fourier 
transformed hopping reads 



V, = 



(18) 
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FIG. 2: Model of an Anderson impurity in the middle of 
the graphene hexagon. Electrons from neighboring can hop 
onto the impurity by the hopping matrix elements Vi. For an 
electron from the vicinity of the Brillouin zone corner K phase 
differences of its wave function at neighboring sites belonging 
to sublattice A (blue marked atoms) are given. The sum of 
these phase factors vanishes. 

As the dispersion does, this coupling vanishes linearly 
when approaching the Brillouin zone corners K and K'. 
As a consequence, any possible Kondo resonance due to 
such an impurity will lead to g ^ 1: Eq. (7) results in 
B ^ 0{Ep) for E-p in this case — a much stronger 
enhancement of the q-factor than for the impurity on top 
of carbon or at a bridge site. 

The origin of this effect can be either understood in 
terms of the Cg symmetry as discussed above or in terms 
of destructive quantum interference in graphene lattice: 
For each state in the vicinity of K (K') the phase of the 
wave function at neighboring sites of the impurity be- 
longing to one sublattice circulates either clockwise or 
counterclockwise around the impurity, as illustrated in 
Fig. 2. The phase of the wave function of atoms in sub- 
lattice A having one common nearest neighbor in sublat- 
tice B circulates around this atom in sublattice B leading 
to £,{k) — > for fc — > X or fc ^ K'. This cancellation 



can also be viewed as a result of Berry phase associated 
with the Dirac points in pristine graphene. The relation 
between vanishing of ^(fc) and topological properties of 
honeycomb lattice was discussed in Ref. 25. Thus, the 
imaginary part ImG{E) vanishes linearly at the conical 
point which leads, taking into account analytical proper- 
ties of G{E) to logarithmic divergence of the q-factor. 

In conclusion, we have addressed the Kondo effect 
in graphene for a realistic d-electron case of Co atom. 
The crystal field of graphene honeycomb lattice splits 
the d-orbitals into two doublets, £'1^2 and one singlet 
Ai states. We have performed the first principles cal- 
culations and found that Ei doublet states are responsi- 
ble for the Kondo effect and for unusual Fano q-factors 
seen in the experiments^^. For the impurity placed in 
the middle of hexagon we have found that the same de- 
structive interference that lead to linearly vanishing DOS 
and to Berry phase is responsible for the anomalously 
large factors in Fano resonance analysis. We thus con- 
clude that nontrivial properties of the Kondo effect in 
graphene are related to masless Dirac fermion spectrum. 
The quantum interference effects, like the Fano effect 
considered here, are extremely sensitive to atomistic de- 
tails such es specific impurity positions, which can be 
checked experimentally. Upon completion of this work 
we learned about the recent preprint by H.-B. Zhuang 
et. al. http://arxiv.org/abs/0905.4548 that addresses 
related questions for the case of a s-wave magnetic im- 
purity. 
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